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ABSTRACT 

Han fc Gouldl 



We present microlensing calculations for a Galactic model based on 
(2003), which is empirically normalised by star counts. We find good agree- 
ment between this model and data recently published by the MACHO and 
OGLE collaborations for the optical depth in various Galactic fields, and the 
trends thereof with Galactic longitude I and latitude b. We produce maps 
of optical depth and, by adopting simple kinematic models, of average event 
time-scales for microlensing towards the Galactic bulge. We also find that 
our model predictions are in reasonable agreement with the OGLE data for 
the expected time-scale distribution. We show that the fractions of events 
with very long and short time-scales due to a lens of mass M are weighted 
by M 2 n(M) dM and M~ l n(M) dM respectively, independent of the density 
and kinematics of the lenses. 

Key words: gravitational lensing - Galaxy: structure - Galaxy: bulge 



1 INTRODUCTION 



A microlensing event occurs when a luminous object is temporarily magnified by a massive 
body, such as a star or dark matter o bject, passing close to the line of sight and acting as 



a gravitational lens. 



Paczvhski 



(1986) advocated searching for microlensing events towards 
the Large Magellanic Cloud in order to detect dark matter in the Gal actic halo. Soon thre e 
separat e collaborations were conducting syst ematic searches: O GLE ((Udalski et al.l ll992h 



EROS (jAubourg et al.lll993f ). and MACHO (jAlcock et al .1 11 9931. between them observing 
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the Large and Small M agellanic Clouds and the Galactic bulge. Other groups such as MOA 



(.e.g. 



Bond et al 



2001) have since joined the search, and thousands of microlen sing events 



Alcock et al 



2000; 



Wozniak et alLoOl : SumikoO^ 



almost all 



microlensing candidates have also been identi- 



have now been detected (e.g. 
towards the bulge. A much smaller number o 
fied t oward the Large Magellanic Cloud (e.g., 

2005^ . One of the main aims of all these observations is to accurately measure the optical 



Alcock et al.ll2000fi and M31 (e.g 



Novati et al 



depth, r - the probability of seeing a microlensing event at any given instant - which can 
provide much information about the structure and mass distribution of the Galaxy and its 
halo. 



Since the first estimates of r bv iPaczyhskil (|l99lh and lOriestl (|l99lh . predictions based 
on increasingly refined models have consistently and significantly disagreed with measure- 
ments based on increasingly l a rge se ts of observational data. However, there are now signs 
of convergence. iHan fc Gouldl ([20031 ) - hereafter HG03 - used star counts from the Hubble 
Space Telescope (HST) to normalise their Galactic model, predicting r = 1.63 x 10 -6 to- 
wards Baade's window (BW), based on lensing of red clump giants (RCGs). They noted 
reasonable agreement with two recent measurements towards the bulge, also based on 
RCGs, of r = 2.0(2.13) ±0 .4 x 10" 6 an d r = 0.94(1.08) ± 0.30 x 10" 6 , from the MA- 



CHO ([Popowski et al. 



20011 ) and EROS dAfonso et al 



The numbers in parentheses are from table 2 of 



2003J) collab orations, respectively. 



Afonso et al 



(2003), who enabled a bet- 



ter comparison between all bulge optical depth measurements to be made by adjusting the 



values for their offset from BW. Now from 7 years of MACHO survey data, iPopowski et al 



( 20041 ) report r = 2.17±° H x 10" 6 at (I, b) = (1.50 ° -2.68°), whi ch is in excellent agreement 



with recent theoreti cal predictions, inc luding the lHan fc Gouldl result. Most recently, from 



the OGLE-II survey 



Sumi et al 



(2005J) find r = 2.37ig;i * !0~ 6 at (l,b) = (1.16°, -2. 75 c 



which is also consistent with the latest MACHO survey value. 

In this paper we generate Monte Carlo simulations of the Galaxy based on HG03. The 
outline of the paper is as follows. £|21 describes the model and theory, and presents our 
results: In fcS.ll we reproduce the HG03 tbw ; an d then compare our predicted r with the 
recent MACHO and OGLE results in various directions. ^3.21 presents maps of optical depth 
and average event time-scale (duration). These maps can be compared with observations in 
any direction. In ^3.31 we predict the event rate as a function of time-scale and compare this 
to the distribution observed by OGLE. In £ 13.41 we show how at both long and short times 



3 

the time-scale distribution is directly related to the lens mass function. We summarise our 
results in 21 



2 THE MODEL 



2.1 Bulge and disc mass models 



Dwek et al 



(1995!) compared various hypothetical mass density models of the bulge to the 
infrared light density profile seen by the Cosmic Background Explorer (COBE) satellite. We 
use the G2 (barred) model from their table 1, with -R max = 5 kpc. The bar is inclined by 
13.4 ° to the Gala ctic centre line of sight, and the distance to the Galactic centre is set at 8 



kpc. 



Dwek et al. 



used 8.5 kpc, so we adjust their model parameters accordingly. The model 
is then normalised by HST star counts (see the end of E J2.2J) . This independent constraint 
can be used to normalise any bulge model. 



For the disc, we use the local disc density model of IZheng et al.l 1)20011 ). as extended to 
the whole disc by HG03. As the disc model is relatively secure (HG03), it will contribute 
only small uncertainties to predictions of the optical depth, so it is not renormalised as for 
the bulge model. 



2.2 Source and lens populations 



The optical depths reported by 



Popowski et al 



( 20041 ) are base d on lensing of RC Gs in the 



bulge, and HG03 assume only bulge RCG sources in their model. ISumi et al.l f)2005r ) observed 
lensing of red giants and red super giants as well as RCGs. We assume that these different 
types of stars follow the same bar density distribution and are bright enough to be seen 
throughout the bar, which corresponds to the case with 7 = in the following eq. (JHJ)- 

Our lens mass function is generated as in HG03. Their unnormalised bulge mass function 
assumes initial star formation according to 



dN/dM = k(M/M } 



brkJ 



(1) 



where M brk = 0.7 M , a = -2.0 f or M > M brk , and a — —1.3 for M < M brk , consistent 



with observations by 



Zoccali et al 



( 20008 ). However HG03 extended this beyond the latter's 
lower limit of M ~ 0.15 M to a brown dwarf cut-off of M ~ 0.03 M Q . We assume objects 
with masses 0.03-0.08 M and 0.08-1 M become brown dwarfs (BD) and main-sequence 
stars (MS) respectively, 1-8 M stars evolve into 0.6 M white dwarfs (WD), 8-40 M stars 
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become 1.35 M Q neutron stars (NS), and anything more massive forms a 5 M Q black hole 
(BH). 



For MS stars we use the mass-luminosity relation of 



Coxl (|l999l ). and take all other lenses 



to be dark. The m odel is then normalise d by comparing extinction-adjusted MS counts to 



HST star counts (Hol tzman et al 



1998) as described in HG03. The same mass function 



and luminosity relation are also used for the disc. Strictly they should be independently 
estimated, but any uncertainties are small compared to others involved as we find disc stars 
account for only ~ 20 per cent of the total number of stars in BW. 



2.3 Kinematic model 



To calculate the event rate, we must also specify the velocities of the lenses, sources and 
observer. The observer velocity vq is assumed to follow the Galactic rotation, so the two 
velocity components in I and b are given by 

220 kms" 1 , 



vo,i = v , m t = Kins v 0th = 0. 

The lens and source velocities in the I and b directions are given by 

V\ = ^rot + ^rand,l) v b = ^rand,b5 



(2) 



(3) 



where the rotation velocity v ro t and the random velocity i> ran d are from lHan & Gouldl (jl 99,1 ) : 
for the disc v TOt = 220 kms" 1 , and for the bar t> rot is given by projecting f max = 100 kms" 1 
across the line of sight according to 



x 



Vrot 



1 kpc^ 

x 



(R < 1 kpc, solid body rotation), 
(R > 1 kpc, flat rotation), 



(4) 



x 2 + y 2 y/ 2 , and the coordinates (x, y, z) have their origin at the Galactic centre, 



where R 

with the x and z axes pointing towards the Earth and the North Galactic Pole respec- 
tively. The random velocity components f ran d,i and f ra nd,b are assumed to have Gaussian 
distributions. For the disc o ih = (30, 20) kms" 1 , and for the bar we use cr XiyiZ = (110, 



Sumi. Ever fc Wozniak 



2003, and 



82.5, 66.3) kms" 1 as found by Irian fc Gouldl (Il 995) using the tensor virial theorem (see also 



Kuiiken fc Rich! 120021 ) . These values should be altered 



slightly as HG03 used a different normalisation. This may affect our results slightly, but it 
is re-assuring that our results based on such a simple kinematic model appear to agree with 
the data quite well (see 



2.4 Optical depth and event rate 



t in any given direction is an average over the optical depths of all the source stars in 
that direction. The optical depth to a particular star is defined as the probability that it is 
within the Einstein radius (see below) of any foreground lenses. Hence mor e distant stars , 



although fainter and less likely to be detected, have higher optical depths (jStanek 



19951 ) 



HG03 accounted for this with the term 7 in the calculation of observed optical depth: 
4ttG Jo 00 dD s D 2 ^p(D s ) J Ds dD d p(D d )D d (D s - D d )/D s 



(r) 



J~dD a Dt' y p(D ! 



(5) 



where D s and D d are the distances to the source and deflector (lens), and p(D s ) and p(D d ) 
are the source number density and lens mass density. RCGs and other bright stars in the 



bulge can be identified independently of 
presented (in a slightly different form) by 



;heir distance, so 7 



0. Eq. © was originally 



Kiraea fc Paczvhskil (J199J), who also derived an 



expression for the lensing event rate V. We give this here in terms of 7, and account for 
variation in lens mass by bringing the term M -1 / 2 inside the integral: 

AG 1 / 2 



r 



dD s D 2 ^p(D s 



C dD d p(D d )v[D d (D s - D d )/MD, 



1 1/2 



SS°dD a D^p{D B ) 
where v is the lens-source relative transverse velocity, 



(6) 



(7) 



and its components in the Galactic I and b coordinates, v\ and t>b, are related to the observer, 
lens and source velocities by 



b = [iyxy - v ) + (v - v s)^j 



(8) 



s / l,b 



where vb and v$ are the deflector (lens) and source transverse velocities; their components 
in the / and b directions are given in eq. (jBJ). 

The time-scale of an even t ty, is defined a s the time taken for a source to cross the 



Einstein radius of the lens ((Paczvhski 



1996): 



V 



UGMD d {D a -D d y 
{ c 2 D s 



1/2 



(9) 
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3 RESULTS 



3.1 Optical depth in MACHO and OGLE fields 

HG03 calculated r = (0.98,0.65,1.63) x 10~ 6 towards BW for bulge, disc, and all lenses 
respectively. Our equivalent values are (1.06, 0.65, 1.71) xl0~ 6 . HG03 noted that the value 
of 7 makes little difference to r for disc lenses, but for bulge lenses r becomes 0.86 x 10~ 6 
when 7=1. We find r = 0.92 x 10~ 6 in this case. Our results for bulge lenses differ by 
7-8 per cent from HG03's due to a slight difference in implementation of the bulge model 
normalisation. We find that allowing MS disc lenses to also act as sources themselves makes 

a negligible difference to the tot al value of r. 

The MACHO measurement (IPopowski et alJhooi of r = 2.17±g;|| x 10~ 6 at (I, b) = 



(1.50°, —2.68°), was obtained from a sub-sample of their observed fields, the 'Central Galactic 
Region' (CGR), which covers 4.5 deg 2 and contains 42 of the 62 RCG microlensing events 
seen. The coordinates (1.50°, —2.68°) are a weighted average position of these fields; the 
unweighted average is (l,b) = (1.55°, —2.82°). Optical depths were also given for a region 
'CGR+3' that contains 3 additional fields, and for all 62 events. In Tabled we compare our 
expected values to each of these results, and to r reported for each of the individual CGR 
fields. 

OGLE's measurement ( Sumi et all2005 ) of r = 2.37±[J;i x KT 6 at (I, b) = (1.16°, -2.75°) 
made use of 32 RCG events, in 20 of their 49 fields, where (I, b) = (1.16°, -2.75°) is the 
weighted average field position, r was also given for each field; we compare our values to all 
of these results in Table El 

Note that any significant disagreement occurs only in individual fields, and that in only 
1 of the 6 fields (MACHO and OGLE) with > 4 events (OGLE #30) does our value lie far 
outside the stated la uncertainty. 

Table El shows the percentage contributions to the total optical depth and event rate 
from the different types of lenses. The disc lenses contribute about 37 per cent of the optical 
depth and a slightly smaller fraction (31 per cent) of the event rate. We see that 62 per cent 
of all events have luminous (MS) lenses, the other 38 per cent are dark (BD, WD, NS and 
BH). The NSs and BHs contribute about 9 per cent of the optical depth but only 4 per cent 
of the event rate. This is because the events caused by stellar remnants on average have 



longer time-scales, and thus they occur 
In their figs. 12 and 14 respectively, 



ess frequent ly. 



Sumi et al. 



(2005J) and IPopowski et all (J2004J) plot 
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Region/field 7V events * (I, b) (°) tmacho(x10- 6 ) r modcl (xl0- 6 ) 



CGR+ 


42 


(1 


50, 


-2 


68) 


9 1 7+0-47 
-0.38 


2 


43 


CGRt 


42 


(1 


55, 


-2 


82) 




2 


33 




■. ) . > 


( j 




_2 


/.->; 


9 07+O.47 
-0.39 


2 




All events 


62 


(3 


18, 


-i 


30) 


1 91+O.21 


1 


32 


108 


6 


(2 


30, 


-2 


65) 


2.04 ±0.92 


2 


31 


109 


2 


(2 


45, 


-3 


20) 


0.58 ±0.41 


1 


96 


113 


3 


(1 


63, 


-2 


78) 


0.55 ±0.35 


2 


31 


114 


3 


(1 


81, 


-3 


50) 


1.19 ±0.74 


1 


87 


118 


7 


(0 


83, 


-3 


07) 


2.85 ± 1.35 


2 


25 


119 





(1 


07, 


-3 


83) 




1 


74 


401 


7 


(2 


02, 


-1 


93) 


5.13 ±2.16 


2 


85 


402 


10 


(1 


27, 


-2 


09) 


3.95 ± 1.50 


2 


89 


403 


4 


(0 


55, 


-2 


32) 


1.16 ±0.66 


2 


83 



Table 1. Comparison of model and MACHO optical depths for the Central Galactic Region (CGR) and individual fields. 
'Number of events seen by MACHO. ^Weighted average (I, b). % Unweighted average (I, b). 



Region/field 


N events 


(l,b 


(°) 


togle(x10- 6 ) 


T-modcl(xl0 6 ) 


All fieldst 


32 


(1.16, 


-2.75) 


ry Qry-\-Q.53 

z - J '-0.43 


2.43 


1 





(1.08, 


-3.62) 




1.87 


2 


1 


(2.23, 


-3.46) 


2.31 ±2.31 


1.85 


3 


4 


(0.11, 


-1.93) 


3.99 ±2.07 


3.20 


4 


5 


(0.43, 


-2.01) 


2.93 ±1.39 


3.09 


20 


1 


(1.68, 


-2.47) 


1.15 ±1.15 


2.54 


21 





(1.80, 


-2.66) 




2.39 


22 


1 


(-0.26, 


-2.95) 


0.79 ±0.79 


2.42 


23 





(-0.50, 


-3.36) 




2.13 


30 


6 


(1.94, 


-2.84) 


8.88 ±3.89 


2.26 


31 


1 


(2.23, 


-2.94) 


2.10 ±2.10 


2.15 


32 


1 


(2.34, 


-3.14) 


0.87 ±0.87 


2.02 


33 


2 


(2.35, 


-3.66) 


9.69 ± 7.38 


1.73 


34 


2 


(1.35, 


-2.40) 


3.80 ±2.69 


2.65 


35 


2 


(3.05, 


-3.00) 


2.99 ±2.19 


1.98 


36 





(3.16, 


-3.20) 




1.85 


37 


2 


(0.00, 


-1.74) 


2.06 ±1.65 


3.39 


38 


2 


(0.97, 


-3.42) 


2.68 ±2.09 


2.01 


39 


3 


(0.53, 


-2.21) 


1.51 ±0.90 


2.92 


15 





(0.98, 


-3.94) 




1.68 


46 





(1.09, 


-4.14) 




1.56 



Table 2. Comparison of model and OGLE optical depths. *Number of events seen by OGLE. ^Weighted average (I, b). 

average optical depths in latitude and longitude strips. We produce similar plots in Fig. [TJ 
with the OGLE and MACHO data points shown. In both sets of strips the model is in good 
agreement with both sets of data. The single data point at negative I is based on only one 
microlensing event, so the discrepancy has low statistical significance. 









Location/type 


of lens 






Bar 


Disc 


BD 


MS 


WD NS 


BH 


Optical depth 


63 


37 


7 


62 


22 6 


3 


Event rate 


69 


31 


17 


62 


17 3 


1 



Table 3. Percentage contributions, to the total predicted r and T, from different types of lens. 
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Figure 1. Average optical depth in latitude (left panel) and longitude (right panel) strips, for —5.5° < b < 5.5° and —5.5° < 
I < 5.5° respectively. The soli d line shows the model pred iction, whi le the open and solid circles are data points from MACHO 
(fig. 14. iPopowski et al.ll2004) and OGLE (fig. 12. ISumi et"aTll2005l) respectively. 



3.2 Maps of optical depth and average event time-scale 

Figs. El and |3] are maps of expected optical depth and average event time-scale. We can 
clearly see higher optical depths and longer time-scales at negative galactic longitude. This 
is due to the inclination of the bar to the line of sight. At positive longitude the bar is 
closer to us, and the line of sight cuts through the bar at a steeper angle. Hence there are 
fewer potential lenses, in either the disc or the bar, between us and any bar source, and so 
r is smaller. Also, objects rotating around the Galactic centre have a smaller component 
of their velocity along the line of sight, so average transverse velocities will be greater, and 
average time-scales shorter. At negative longitude, the line of sight passes through more 
of the disc and cuts the bar at a shallower angle. Hence we see higher optical depths and 
smaller transverse velocities, and thus longe r average time-scales . 



We wish to compare our maps to others. 



Evans &: Belokurovl ((20021 ) produced red clump 



optical depth maps for three Galactic models, but while two of these appear similar to ours, 
they do not agree. One of those models was also used to make a time-scale map, which is 
quite different to ours. This is not surprising since, as well as using a different mass model, 
their mass function, velocities and velocity dispersions were also different. (In fact their 
timescale map has two sets of contours, to show the effects of including and excluding bar 
streaming. Without streaming their mean timescales are much shorter than ours, and with 
it they are greater by a factor ~ 3, much longer than ours. Such a large variation is puzzling, 
and we are cautious about comparing th e ir ma p to ours). 

In their fig. 16, iBissantz fc Gerhard! ([20021 ) presented an optical depth map for RCG 
sources, with a bar angle of 20°. For b < — 3° it appears quite similar to ours, but moving 
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-10 



Galactic longitude I (°) 

Figure 2. Map of expected optical depth. The MACHO and OGLE fields are shown by the large and small grey boxes, 
respectively. For the MACHO fields, the Crosshatch pattern indicates the CGR subset listed in Table Q For OGLE, 
the Crosshatch pattern denotes those fields not listed in Table [2] The small square indicates BW. Contour levels are at 
(0.1, 0.2, 0.3, 0.5, 1, 2, 3, 4, 5) x 10" 6 . 




Galactic longitude I (°) 



Figure 3. Map of expected average event time-scale. The MACHO and OGLE fields are shown by the large and small grey 
boxes, respectively. For the MACHO fields, the Crosshatch pattern indicates the CGR subset listed in Table HI For OGLE, the 
Crosshatch pattern denotes those fields not listed in Table l2l The small square indicates BW. Contour levels are at 22.5, 25, 
27.5, 30, 32.5, 35, 37.5, 40 and 42.5 d. 



towards the Galactic centre r climbs far more steeply than in our map. This is best seen by 
comparison to their fig. 17, where they plot r as a function of b, for / = 3.9°. This is s hown in 
Fig. 0] with equivalent profiles for our model. We see how rapidlv lBissantz fc Gerhard! ' s profile 
diverges from ours towards 6 = 0°. We also see that changing the bar angle in our model 
from 13.4° to 20° does not explain this difference. Instead it is probably due to the density 
in their bulge mass model increasing much faster towards the mid-plane. The observational 
data for the mid-plane are limited due to heavy extinction, and so mass models are not 
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30 



2::; 



20 



O 15 



10 - 



Bar angle = 13.4 
Bar angle = 20-0* 
Bissantz & Gerhard (2002) 



-8 -6 -4 -2 

Galactic latitude b (°) 

Figu re 4. r as a function of b, for I = 3.9°. Our model slope is almost identical for bar angles of 13.4° and 20°. The profile of 
iBissantz fc Gerhard! J2002T1 diverges from ours at b fa —3°, increasing rapidly towards the mid-plane (b = 0°). 



well-constrained in this region. Given the difficulty in obtaining any measurement of r at 
small latitude, it is difficult at present to test either profile there. 



3.3 Time-scale distributions 



Fig. El shows the event rate as a function of time-scale towards the OGLE coordinates 
(l,b) = (1.16°, —2.75°), for bar (thin line), disc (dashed line) and all (bold line) lenses. 
There is good agreement with the asymptotic power-law tails dT/d(lo R t E ) oc t E 3 , t E ~ 3 
for very short and long time-scales, respectively (jMao fc Paczvhskilll99fi[ ). The disc lensing 
events have an average time-scale of 26.3 d, slightly longer th an the bulge lensing events' 



Kiraea &: Paczvhskil (J199J). The average time-scale for 



average of 25.7 d, as also found by 
all events is 25.9 d. 

In Fig. El we renormalise our time-scale distribution (for all lenses) and compare it t o 
that seen by OGLE, as corrected for detection efficiency (s ee fig. 14 in 



Sum et al. 



We do not compare to the time-scale distribution seen by Popowski et ah ( 2004^ - they 
assumed that the effect of blending on RCG sources is negligible, but ISumi et al.l (2005) 
found ~ 38 per cent of OGLE-II events with apparent RCG sources were really due to 
faint stars blended with a bright companion. Fortunately, they also showed that blending 



200,4 



has little effe ct on estimates o : 
also made by 



r du e to partial cancellation of its different effects, a point 



Popowski et al.l (j2004l ). However, time-scale distributions will be significantly 



shifted towards shorter events. As a result, the MACHO time-scale distribution (not shown) 
has a significant excess at short time-scales compared with our model. 



11 



T 




og t E (days) 

Figure 5. Predicted microlensing event rate as a function of time-scale, towards (l,b) = (1.16° , —2.75°). The bold line represents 
all lenses. The thin and dashed lines represent the bar and disc lenses, respectively. The two dotted lines are asymptotic tails 
dr/d(log t E ) oc t E s , t E - 3 for very short and long time-scales, respectively. 



6 r 

5 r 




t E (days) 

Figure 6. Microlensing event rate as a function of time-scale, towards (I, b) = (1.16°, —2.75°). The solid line shows the model 
prediction, and the OGLE observed distribution (corrected for detection efficiency) is shown as a histogram. 

Our time-scale distribution shows reasonable agreement with OGLE's. The Kolmogorov- 
Smirnov (KS) test shows that the predicted and observed distributions are consistent at a 
w 52 per cent confidence level. Our average time-scale of 25.9 d is in excellent agreement 
with OGLE's corrected average of 28.1 ±4.3 d. Our median and quartiles are (19.2, 11.2, 
31.7) d, respectively. 

The event time-scale distribution from the data still has large uncertainties due to the 
limited number of events. It is apparent that the data have not yet reached the predicted 
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asympototic behaviour at short and long time scales, so a more stringent test on the model 
is not yet possible. 



Bissantz. Debattista fc Gerhard! ([2004L see also 



Peak 



1998) have also modelled the time- 



scale distribution. They reproduced that from MACHO's 99 DIA events ([Alcock et al 



20001) 



centred at (l,b) = (2.68°, —3.35°). However, both distributions are clearly shifted towards 
short time-scales compared to our model prediction in the same direction 1 (this is not shown, 



as it is very close to the solid line in 



systematics of blending (jSumi et al 



FigJ 



200E 



. Although the DIA method is less prone to the 
it is still possible that the MACHO DIA time- 



scal e distribution i s somewhat affected. The most important difference between our model 



and 



Bissantz et al. 



's is that in order to match the data at short time-scales, they adopted a 
Schecter mass function, n(M) oc M~ 2 35 for M < 0.35 M down to 0.04 M Q , steeper than our 
mass function, n(M) oc M~ L3 for M < 0.7 M Q . As a result, their median lens mass is much 
smaller than ours (0.11 M vs. 0.35 M , weighted by event rate). The different kinematics 
may also have a noticeable effect on the timescales, but their more realistic dynamical model 
does not allow a simple comparison to be made. 



3.4 Fractional contributions to event rate — mass weightings 

Fig. [7| shows the fractional contributions to the total event rate, as a function of event 
time-scale, for the different types of lens (BD, MS, WD, NS and BH) as indicated. At 
short time-scales (t-E<4 d), the brown dwarfs dominate the event rate, while at long time- 
scales (£e ^ 100 d), the stellar remnants become increasingly important. There is asymptotic 
behaviour at both long and short time-scales. We find that the fractional contribution from 
a lens of mass M is weighted by M 2 n(M) dM and M _1 n(M) dM, respectively. In the 
Appendix we de rive these weighti ngs from eq. ©• (The scaling at long event tails has already 
been derived bv lAgol et al.ll2002l ). Tabled shows that direct calculation of these asymptotic 
fractions from the mass function gives results that clearly agree with the trends in Fig. [7| 
These weightings are independent of the density and kinematics of the lens population, and 
hence provide valuable information about the lens mass function. 



1 At first glance all three distributions may appear to be similar. However, whereas we define the event timescale as the 
Einstein-radius crossing time (see 112.41 . MACHO plot the diameter-crossing time, a factor of 2 difference. 
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Time-scale 


BD 


MS 


WD 


NS 


BH 


Long 


0.53 


44 


20 


12 


24 


Short 


72 


27 


1.5 


0.078 


0.0032 



Table 4. Percentage contributions to the total predicted event rate, at long and short time-scales, from the different types of 
lens. 
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Figure 7. Fractional contributions to total expected event rate, as a function of event time-scale, from BD, MS, WD, NS and 
BH lenses as indicated. The solid and dashed lines represent the bar and disc lenses, respectively. The asymptotic fractions at 
long and short time-scales are a function of the lens mass only (see text). 

4 SUMMARY 

In th is paper, we have used a simple Galaxy model normalised by star counts ( Han &: Gouldl 
20031 1 to predict the microlensing optical depth. Combined with simple kinematic models, we 



also predict maps and distributions of the time-scale distributions. We have shown that the 
fraction of long and short events contributed by a lens of mass M is weighted by M 2 n(M) dM 
and M _1 n(M) dM respectively. If the tails of this distribution can be accurately determined 
from observations, we have a direct probe of the lens mass function. 

It is remarkab le that this emprically-normalised model based on the COBE G2 model 



fDw ek et al 



1995) sho ws good agreemen t wi th data recently publi shed by the MACHO and 



OGLE collaborations ( Sumi et al 



2005 and 



Popowski et al. 



p. 

2( 



2004J) for the optical depth in 



various Galactic fields, and its trends with I and b. Our maps of optical depth and average 
event time-scale cover a large area of the sky, and can be compared to future determinations 
of r in similar areas when they become available. The expected distribution of the event 



time-scale a 


so ap 


(Sumi et al. 


2005|) 



However, the numbers of microlensing events used (42 and 32) in the 



recent MACHO and OGLE analyses are still small, so the test on the models is not yet 
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stringent. When the much larger database of microlensing events (~ thousands) is analysed, 
then a full comparison with the models will become much more discriminating. 
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APPENDIX A: EVENT RATE WEIGHTINGS AT LONG AND SHORT 
TIME-SCALES 

As described in H'SA\ the microlensing event rate shows asymptotic behaviour at both long 
and short time-scales. We show here that this is directly related to the lens mass function, 
specifically, the fractional contributions are weig hted by M 2 n{M) dM and M -1 n(M) dM, 
at very long and short time-scales respectively. 

The event rate is given by eq. (jHI) • However, as the mass dependence of the asympototic 
behaviour is the same for sources at different distances, we shall ignore the source distance 
dependences here. Therefore for a source at distance D s and a population of lenses each with 
identical velocity v and mass M, the event rate is given by 

4^1/2^ ( D d (D s -D d ) \ 1/2 ^ 

T = —J p{Dd)v { MD S ) dDd ' (A1) 

where p(-Dd) is the lens mass density at D d . 
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In reality, v and M both vary. The velocity probability distribution, p(v) dv, can usually 
be approximated by a two-dimensional Maxwellian distribution 

p(u) dv = — exp(-v 2 /2a 2 ) v dv, (A2) 



(T 



where a is the velocity dispersion. For constant M, the factor p(D d ) in eq. (jAl|) is simply 
the total mass density. When M varies, the event rate depends on the lens mass function, 
i.e. on how the total mass is partitioned into lenses of different masses. We assume that this 
is the same everywhere. The mass density for lenses with M — ► M + dM can be written as 
a product of i(D d ) n(M) dM, where i(D d ) indicates the distance dependence and n(M)dM 
is the number density of lenses between M — > M + dM. 

Integrating over the mass and velocity distributions and using the fact that p(M)dM = 
Mi(D d ) n(M) dM, we obtain 

r = 2A 1/2 J°* D cS 1/2 t(D d ) dD d J n(M)M l/2 dM J v p(«) dv, (A3) 

where A = 4G/c 2 and D eS = D d (D s — D d )/D s . We can now rewrite the time-scale equation 
(eq.EJ) 

r E [AGMD d (D s -D d )\ 1/2 



v 



-i 



V c 2 D s 
A l / 2 M l / 2 D cS l/2 



v 

The typical transverse velocity is ~ a, and this defines a characteristic time-scale as 

A l/2 M l/ 2jDeff l/2 
t a — . 

a 

The short and long tails satisfy t E t a and t E 3> t a , respectively. 



(A4) 



(A5) 



Al Behaviour at long time-scales 

As can be seen from eq. (|A4|) . the long time-scale events occur when the lens and source 
both move approximately parallel to each other and perpendicular to the line of sight. In 
this case, the transverse velocity is close to zero (v a) and the time-scale becomes long. 
For events with time-scales longer than ti ong (^ > t a ), the transverse velocity must satisfy 

A 1/2 M 1/2 D 1/2 

v < — < a. (A6) 

^long 

The exponential factor exp(— v 2 /2a 2 ) approaches unity, and so we have 
r(>t long ) = —I D eS 1 /H(D d )dD d J n{M)M 1 ' 2 dM jf tlong v 2 dv 
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2A 2 r D s r 
= q 2+ 3 / D cS H(D d )dD d n(M)M 2 dM. (A7) 

JO" *long JO J 

Therefore, for long time-scale events, the event rate follows a power-law as a function of 
time-scale, wit h a normali s ation that depends on the mass function, oc M 2 n(M) dM, as 



also derived by 



Agol et al 



(2IM) 



A2 Behaviour at short time-scales 

Re-expressing equation (|A4|) in terms of x = D d /D s , we have 
/4-OM \ 1 I 2 

t E =(^-x(l-x)D B ) v-\ (A8) 

Very short events occur when the lens is very close to either the source or the observer, i.e., 
when x — > 1 or x — > 0. The asymptotic behaviour is the same for 2; < 1 and 1 — x <C 1, 
so we concentrate here on the case when i<1, x(l — x) ~ x. So for events shorter than a 
given time-scale t s hort(^ to), we must have 

x < < A9 > 

Equation ()A3|) can then be re- written in terms of x: 

r(<t 9 hort) = 2A 1 ' 2 j° S [x(l-x)D s ] l l 2 i(D A )dD d J n{M)M l ' 2 dM J vp(v)dv. (A10) 

Changing the integration variable to x, and with i < 1, f(xD s ) ~ f(0), we obtain for the 
first integral 



T(<t 



short ) = 2A 1 ' 2 J AMDs x l/ H{xD s )D^ 2 dx J n{M)M 1 ' 2 dM J vp{v)dv 

1 r: '-l(0) / M~ l n(M) dM [ v 4 p(v)dv. (All) 
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Therefore for short time-scale events, the event rate follows a power-law as a function of the 
time-scale, with a normalisation that depends on the mass function, oc M -1 n(M) dM. 



